Shape Selection in Chiral Self- Assembly 
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Many biological and synthetic materials self- assemble into helical or twisted aggregates. The 
shape is determined by a complex interplay between elastic forces and the orientation and chirality 
of the constituent molecules. We study this interplay through Monte Carlo simulations, with an 
accelerated algorithm motivated by the growth of an aggregate out of solution. The simulations 
show that the curvature changes smoothly from cylindrical to saddle-like as the elastic moduli are 
varied. Remarkably, aggregates of either handedness form from molecules of a single handedness, 
depending on the molecular orientation. 



A wide range of materials self-assemble into chi- 
ral aggregates. These include biological materials and 
their synthetic analogues, such as carbohydrate- [j| and 
amino-acid-based amphiphiles 2], peptides 0, yjj, di- 
acetylenic lipids 0, p , gemini surfactants Q , and mul- 
ticomponent mixtures in bile @, Ej . Aggregates form as 
bilayer membranes in solution, and grow into a range 
of chiral shapes, including tubules with "barber-pole" 
markings, helical ribbons with cylindrical curvature, and 
twisted ribbons with Gaussian saddle-like curvature [HI . 
They form through a variety of pathways, including reor- 
ganization of large spherical vesicles and growth directly 
out of solution pj],[T3. In the latter process, aggregates 
twist in a helix while they increase in length. 

While one might expect that macroscopic handedness 
of an aggregate is tied to the molecular-scale chirality of 
its constituents, experiments have shown that the rela- 
tionship is not always that simple. In the early stages 
of self-assembly, diacetylenic lipids can form both right- 
and left-handed helical ribbons, even though the equi- 
librium state has tubules of a uniform handedness [l2(. 
In a biological analogy, chirality reversal at the organism 
level arises occasionally via genetic mutations, and has 
been observed in plants, molluscs, birds, and mammals, 
even though the proteins and other molecules within each 
organism retain their normal chirality [IH . H3 . HEj . 

Experiments have also shown that the relationship be- 
tween material properties and the macroscopic shape of 
chiral aggregates is surprisingly complex. For example, 
charged gemini surfactants with chiral counterions ex- 
hibit a transition between helical ribbons with cylindri- 
cal curvature and twisted ribbons with Gaussian curva- 
ture as a function of molecular chain length Simi- 
larly, mixed bilayers of saturated and diacetylenic phos- 
pholipids show transitions from micron-scale cylindrical 
tubules to twisted ribbons, and then to nanometer-scale 
tubules, as a function of temperature |l6j . 

To understand the shape of chiral aggregates, and to 
design new structures for nanotechnology, we must learn 
how to control aggregate size and shape by adjusting the 
chemical composition and the conditions under which 
self-assembly occurs. An important step toward this 



goal is to understand how the membrane's elastic moduli, 
chirality, and tilt orientation together control the struc- 
ture of the resulting aggregate. Earlier theoretical mod- 
els have studied this by using continuum elastic theory 
to calculate the energy of twisted membranes in simple 
geometries pjj. Such models show that chiral interac- 
tions, coupled with molecular tilt, lead to the formation 
of cylindrical tubules or helical or twisted ribbons. The 
limitation of this analytical approach is that it considers 
only idealized structures, and does not explore the full 
range of potential shapes or the mechanism of shape tran- 
sitions. Finite-element modeling of membranes and vesi- 
cles would allow free exploration of shapes, but no such 
models have taken the membrane's molecular tilt and chi- 
ral interactions explicitly into account. Molecular-scale 
simulation of tilted bilayer membranes is another possi- 
ble approach , but available system sizes are too small 
to allow studies of shape selection. 

In this paper, we present a new mesoscale simulation 
approach to explore the roles of elasticity, chirality, and 
molecular tilt in chiral self-assembly. We develop a simu- 
lation technique motivated by the growth of an aggregate 
out of solution, which allows the aggregate to select its 
shape through a Monte Carlo process. This approach 
provides an explicit visualization of a range of possi- 
ble shapes, and shows that the aggregate shape varies 
smoothly from cylindrical to Gaussian curvature as the 
elastic properties of the membrane are varied. Moreover, 
it shows that the handedness of the aggregate depends on 
the orientation of molecular tilt as well as on the handed- 
ness of the chiral interactions, demonstrating one mech- 
anism by which both right- and left-handed structures 
can form from a material of uniform molecular chirality. 
Thus, our simulations may explain a wide range of ex- 
perimental findings: the formation of helical and twisted 
ribbons, the transitions between these shapes in gemini 
surfactants Q and lipid mixtures 10], and the observa- 
tion of chiral aggregates with both helical senses in the 
same lipid [l2|. 

In our mesoscale approach, a chiral elastic ribbon is 
modeled as a tethered membrane composed of interac- 
tion sites connected by elastic springs to form a trian- 
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FIG. 1: (Color online) Model of a chiral elastic ribbon as a 
tethered membrane composed of interaction sites connected 
by elastic springs to form a triangular lattice. The local 
molecular orientation at each site is described by a director, 
as indicated by the arrows. 



gular lattice, as shown in Fig. 1. Each interaction site 
represents a small "patch" of membrane. The membrane 
can deform in three-dimensional (3D) space, so that each 
site has position r^, and local upward normal n^. The lo- 
cal molecular orientation at each site is described by the 
director d^, a vector of unit length. This approach is 
analogous to earlier studies of tethered membranes pjj? 
but here the membrane has the additional features of chi- 
rality and molecular orientation. It is also analogous to 
mesoscale models of smectic liquid crystals, which have 
chirality as well as a molecular orientation defined on a 
lattice |20j . but here the membrane is free to change its 
shape in 3D. In our initial exploration of the model, we 
constrain the director on each site at a typical experi- 
mental tilt angle of 30° with respect to the local upward 
normal n^, with the tilt oriented along the ribbon's long 
edge, which is also along a bond between neighboring 
sites. In future work we will remove this constraint and 
consider the possibility of modulated membrane tilt. 

To investigate shape selection, we must define a model 
for the elastic energy of a chiral membrane that takes 
into account the microscopic chirality of its constituent 
material. We use a lattice expression analogous to the 
continuum energy for chiral membranes investigated in 
earlier work [Tjj. It includes a chiral elastic term as well 
as bend, stretch, and hard-core contributions: 
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In the first term of this expression, the chiral elastic pa- 
rameter A gives an energetic preference for twisting the 
directors and dj away from parallel alignment. Re- 
versing the sign of A reverses the favored chiral twist of 
the director field. In the second term, the parameter 
^bend controls the energy of bending the membrane so 
that the local normal vectors on sites i and j are not par- 




FIG. 2: (Color online) Three views of a helical ribbon, show- 
ing the 3D structure. The ribbon has grown to a size of 
16 x 216 sites. 



allel. In the third term, the Hooke's law constant if s tretch 
gives the energy of separating sites i and j away from 
their equilibrium separation Rq. In those three terms, 
the sum on (i, j) refers to all nearest-neighbor pairs of 
sites in the membrane. The final term in the energy rep- 
resents a hard-core potential between any pair of sites, 
with Uhc(r) = for r > a, and ?7hc(0 = oo for r < a, 
where a = |i?o is the hard-core diameter. That term 
prevents the membrane from intersecting itself. 

For Monte Carlo simulations of the membrane, we need 
a simulation algorithm that will come to equilibrium in 
a reasonable time. For a first attempt, we initialized the 
system as a long, flat ribbon, and allowed it to evolve via 
off-lattice single-particle moves with the Metropolis algo- 
rithm. However, we found that this process evolves much 
too slowly to reach a twisted shape. Because the twist 
must start at the ribbon's narrow ends and diffuse inward 
toward the middle, there is a barrier against twisting a 
ribbon, and the equilibration time grows sharply with the 
ribbon length. 

To develop an alternative simulation approach, we are 
motivated by one of the experimental pathways for the 
formation of chiral aggregates. When ribbons grow out of 
solution, they wrap into a helix as they grow, and hence 
they have no barrier against twisting. For that reason, 
we simulate a growth process, with particles attaching 
onto one edge of the model aggregate. The concept of 
simulating a growth process has found wide application 
in the study of polymers [2l|, [22, but this is its first 
application to chiral membranes. While this simulation 
approach mimics an observed growth mechanism, it is 
essentially a mesoscopic Monte Carlo algorithm, which 
does not take into account any of the detailed kinetics 
associated with molecular aggregation. 

We start the system as a small ribbon of size 16 x 16 
sites, and then add new rows one at a time up to a max- 
imum size of 16 x 216 sites. After each row addition, the 
last 8 rows are equilibrated for 10,000 Monte Carlo steps 
per site, and then the whole ribbon is equilibrated for an 
additional 5000 Monte Carlo steps per site. The temper- 
ature is kept near zero so that the Monte Carlo process 
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FIG. 3: (Color online) (a) Series of ribbons grown with the 
parameter if stretch = 2, 4, 6, 8, and 10, showing the crossover 
from twisted ribbons with Gaussian saddle-like curvature to 
helical ribbons with cylindrical curvature, (b) Plot of the 
Gaussian curvature for the five ribbons, calculated at the cen- 
tral point of the ribbons, in units of the hard-core diameter 
a -2 . This plot confirms that all of the structures have some 
Gaussian curvature, and that the level of Gaussian curvature 
changes gradually as a function of the elastic parameters. 



is essentially an energy minimization, which allows com- 
parison with previous analytical results. In future work 
we will address the role of thermal fluctuations in mem- 
brane shape. An example of a chiral membrane grown 
in this fashion is shown in Fig. 2 . an d animations of the 
growth are presented in EPAPS 231. In the fi smres and 
animations, colors are a guide to the eye to visualize the 
3D structure. 

We use this simulation approach to investigate how 
the final equilibrium shape of the aggregate depends on 
membrane properties. Of the three elastic parameters 
^stretch, ^bend, and A, it is necessary to vary two since 
only their ratios are significant at low temperature. We 
begin by examining the effect of varying the stretch mod- 
ulus if stretch over a range from 2 to 10 in arbitrary units, 
while ifbend = 1 and A = 1 are held fixed. We ob- 
serve the gradual structural transition shown in Fig. 3(a). 
For low if stretch, we find a twisted ribbon with Gaussian 
saddle-like curvature, while for high if stretch, the shape is 
a helical ribbon with cylindrical curvature. This transi- 
tion is physically reasonable, because Gaussian curvature 
requires stretching the membrane, which is suppressed 
by if stretch- The small Gaussian curvature that remains 
in helical ribbons for high if s tretch resembles the ripple 
structure that has been predicted theoretically J24|. 

We have done an analogous study of the effects of vary- 
ing the bend modulus if bend- We fix if stretch = 6 and 
A = 1, and vary if bend over the range from 1 to 2.4. 
The simulation results (not shown) give a similar gradual 
transition from helical ribbons at low if bend to twisted 



FIG. 4: (Color online) Two ribbons grown with the same elas- 
tic parameters, but with different orientations of the director 
representing the molecular tilt. In the left case, the tilt is 
aligned along the long axis of the ribbon, but in the right 
case it is aligned along the short axis. This change in the tilt 
direction reverses the handedness of the ribbon, and greatly 
changes the radius and pitch. 



ribbons at high if bend- 

Each structure is characterized by a radius, a pitch, 
and two curvature parameters. Through the structural 
transition shown in Fig. 3(a), all of these measures show 
a continuous and gradual transition. In particular, note 
that none of the structures are ideal twisted or heli- 
cal ribbons with pure Gaussian or cylindrical curva- 
ture. Rather, each structure involves a combination of 
Gaussian and cylindrical curvature, and the combination 
changes gradually as a function of if stretch- We can cal- 
culate the Gaussian curvature at any interaction site by 
adding up the six incident bond angles 0i at that site. If 
the Gaussian curvature is zero, then the six bond angles 
add up to 27r. Hence, the deviation of the sum away from 
2tt is a simple measure of the Gaussian curvature |25ll2q : 
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where a$ are the areas of the triangles adjacent to the 
site. Figure 3(b) plots the Gaussian curvature of each 
ribbon shown in Fig. 3(a), evaluated at the central point 
of the ribbon, as a function of if stretch- The variation 
in Gaussian curvature is very substantial, ranging over 
a factor of 5. However, it does not occur as a sharp 
structural transition, but rather as a smooth crossover. 

To further investigate this model, we explore the con- 
sequences of changing the molecular tilt direction in the 
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membrane. In the original simulations, the director was 
constrained to point along the ribbon's long edge, along 
a bond between neighboring interaction sites (as in a 
smectic-J liquid crystal). We now rotate it by 90° to 
point perpendicular to the long edge, halfway between 
two nearest-neighbor bonds (as in a smectic-F liquid 
crystal). Although the sign and magnitude of the mem- 
brane's chiral interaction parameter A remain unchanged, 
the handedness of the aggregate switches from right- to 
left-handed, as shown in Fig. 4. The pitch and diameter 
of the aggregate also change greatly. 

While somewhat counter-intuitive, the handedness re- 
versal is physically reasonable. Earlier theoretical work 
showed that the curvature direction favored by molecu- 
lar chirality is 45° from the tilt direction [Tjj, so rotat- 
ing the tilt direction by 90° should change the curvature 
direction by 90°, giving a handedness reversal. (As a 
thought experiment, we can imagine cutting a cylindri- 
cal tubule to form a helical ribbon. If we cut it along the 
tilt direction, the ribbon has one handedness. If we cut 
it perpendicular to the tilt direction, the ribbon has the 
opposite handedness.) This result shows that aggregates 
may form with either handedness even from enantiomer- 
ically pure material, simply by variation of the molecular 
tilt direction with respect to the edges. This could be 
a possible explanation of the experiments on conflicting 
handedness mentioned earlier |l2|. 

The change in pitch and diameter is a more subtle 
effect, which we cannot attribute to the change in tilt 
direction relative to the ribbon's long edge. Rather, it 



is probably due to changing the tilt direction relative 
to the nearest-neighbor bonds in the membrane. This 
change can greatly modify the elastic properties of the 
membrane, which would alter the geometric parameters 
of the resulting aggregate. This result shows that ag- 
gregates can have different pitch and diameter even if 
they are composed of the same material, just by variation 
of the tilt direction with respect to the crystallographic 
axes. This could be a possible explanation of the different 
geometric structures seen in bile 0, . 

In conclusion, we have developed an accelerated ap- 
proach for studying the shapes of chiral aggregates 
through Monte Carlo simulation of a growth process. 
The simulation results show that a short-range chiral in- 
teraction leads to the formation of complex chiral shapes, 
which undergo a gradual crossover between Gaussian 
saddle-like curvature and cylindrical curvature as a func- 
tion of the interaction parameters. The handedness of 
the macroscopic shapes is determined by the interplay be- 
tween the microscopic chiral energetic parameter and the 
orientation of the molecular tilt. Thus, this work shows 
the range of shapes that occur in experimental studies of 
chiral self-assembly, provides a possible explanation for 
experiments on conflicting handedness, and demonstrates 
an approach for future simulations based on even more 
detailed microscopic models of chiral molecular structure. 
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